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Abstract 

In this paper we consider two-point unilateral Markov fields on a two-dimensional 
lattice as considered by Pickard jT3], Galbraith and Walley [1J[3]. We show that, 
under various ergodicity conditions, they can be perfectly simulated in the station- 
ary state on any finite window. The techniques which are used connect perfect 
simulation with oriented percolation through suitable coupling constructions. 
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1 Introduction 

In this paper we consider a particular class of random fields on the two-dimensional lattice 
N +2 and their stationary extensions to Z 2 , the so-called unilateral random fields. These 
fields appeared in the literature as models for crystal growth [Tj\; later they have been 
proposed also for image analysis applications [H] . Unilateral random fields can be defined 
in general w.r.t. an oriented acyclic graph structure (as in [TU] . where they are called 
Bayesian networks): for the sake of simplicity we will consider the so-called two-point 
models on Z 2 , which corresponds to the nearest neighbor structure where oriented edges 
connect with (i + 1, j) and + 1). However the reader will recognize that our 

basic idea works for more general translation invariant graphs. 

A unilateral two-point field model is constructed through a transition kernel. We 
assume that the state space E is Borel, i.e. a Borel subset of a Polish space endowed 
with its Borel a-algebra. A (unilateral two-point) transition kernel K on E is a mapping 
K : E 2 — > V(E), where V(E) is the space of probability measures on E, with the property 
that for any measurable set B C E the function K(B\x,y) = K(x,y)(B) is measurable 
in the pair (x, y) G E 2 . An E- valued unilateral field X x = {Xfj G E, G N +2 } can be 
constructed from K for any choice of the boundary conditions (b.c.'s) x = {x^q, Xq^, % G 
N + } . The construction is accomplished by specifying consistently the law of X* m n in any 
finite box 

A m>n = {1,2,..., to} x {l,2,...,n} (1) 
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as being 

„ m n 

P(XX m , n GB) = / nil K{dx id \ Xi ^ (2) 
t=i i=i 

where B is any measurable set in E mxn . A natural way to simulate this field is to produce 
the variables Xf^ in any sequential (total) order such that each site comes after its 
parents (i — l,j) and (z, j — A total order with this property will be called increasing in 
the sequel. The reader will immediately notice the similarity with discrete time Markov 
chains. 

Next we randomize the b.c.'s, say with law fi. We define the law of the unilateral field 

as 

P(XL,„ eA)=J P(Xl m n G AMdx). (3) 

We say that the b.c.'s are Markovian when the horizontal boundary X. = {Xi^,i = 
1,2,...} and the vertical boundary Xq,. = {Xq^, i — 1, 2, . . . } are Markov chains which are 
conditionally independent given a common starting value X^q (which can be considered 
as the value of the field at the origin). For Markovian b.c.'s the sequential simulation 
of the field X M is easy, using any increasing order of the sites of N 2 (including boundary 
sites). 

Pickard studied the problem of determining laws fi such that the corresponding 
field X^ 1 is stationary, i.e. invariant under translations in N 2 . Stationary unilateral fields 
can be extended to Z 2 by Kolmogorov's theorem. When E is finite one can prove that 
stationary unilateral fields can always be constructed. Using again Kolmogorov's theorem, 
it suffices to take compatible convergent subsequences of the sequence of averages 

A*w(0 = ^ E P ( X V^-), L = l,2,... (4) 

teA L ,L 

in any fixed finite box A m>n , and use standard compactness arguments. 

For Markovian b.c.'s, Pickard established sufficient conditions for stationarity of a 
unilateral field, which become also necessary in the binary case. Pickard's conditions 
are formulated in terms of the joint distribution of the 4-tuple of random variables 
(Xq o, X^ Q , Xq 1( X^) (see [IH]), hence they involve K and //. On the other hand, Gal- 
braith and Walley have proved that in the binary case any positive two-point unilat- 
eral kernel K has a unique Markovian boundary law jjr = \i under which the columns 
(Xq o, Xi ) and (X^, X^) have the same law, and also the rows (Xq Q , Xq X ) and (Xf , X^) 
have the same law 0j. Since these two conditions are necessary for a stationary field, they 
could also reformulate Pickard's conditions directly in terms of K. 

However not all stationary unilateral fields have Markovian boundary laws, hence they 
are not necessarily easy to simulate sequentially. The goal of this paper is to construct 
simulation algorithms for the class of stationary unilateral fields which is introduced next. 

We say that a unilateral field X x with kernel K is ergodic if there exists a stationary 
field X' = {X^ G E, i,j = 0, 1, 2, . . . } such that for every b.c.'s x and every pair (m, n) G 
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N 2 

Hm \\£(X_f Kl)+Am J - £(X> Am J\\ = (5) 

where we use C(Y) to denote the law of a random vector Y, and || • || is the total 
variation. If the limit in (pj) is uniform in the boundary conditions x we say that the 
kernel is uniformly ergodic. In the binary case Galbraith and Walley |E] gave a quite 
involved sufficient condition for the uniform ergodicity of K. 

We call the field X' appearing in (pj) the equilibrium field of the kernel K. It is easy to 
see that X' is a stationary unilateral field with kernel K, and if it exists it is necessarily 
unique. Since its boundary law // could be very hard to compute (even in the binary 
case), it is not easy to simulate X'a even on a small box A. In fact, by © we know 
that we can only approach £(X'a) by shifting the box sufficiently far away from the 
boundary. The problem is analogous to that of sampling a stationary Markov chain on a 
finite window when the stationary distribution is not available (as it happens in MCMC 
simulations); but in this case, in order to simulate a box, we need to determine the whole 
joint equilibrium law of its boundary. 

Under a rather strong minorization condition on K (Assumption 12. lj) for a general 
Borel state space E, we show in the next two sections how to implement a simulation 
algorithm which produces a sample on any finite box, exactly distributed as the equilib- 
rium field. As a byproduct, we establish uniform ergodicity of K. For one-dimensional 
discrete-time Markov processes the first algorithm of this type was the CFTP algorithm 
of Propp and Wilson ^H] , see also [B] . This work stimulated a wide interest toward what 
is now called perfect simulation, see |12j . Our algorithm is based on an idea introduced 
by Murdoch and Green ^3] for discrete-time Markov processes. It consists in coupling 
the whole family of unilateral fields, with all possible boundary conditions, by using an 
underlying auxiliary Bernoulli field. We show that the dependence from boundary values 
is propagated only along increasing open paths in such a Bernoulli field. In this way 
results from oriented percolation for two-dimensional Bernoulli fields can be used in order 
to ensure that the propagation stops with probability 1 when the boundary moves far 
away from the region to be simulated. In order to simulate a square box of side L, the 
algorithm requires, in addition to the L 2 variables of the box, an average of O(L) addi- 
tional random variables of the field. For other graphical constructions used in the study 
of ergodicity see e.g. j2]. 

In Section 01 in order to relax Assumption 12.11 we study a more general class of 
algorithms working on blocks of adjacent sites lying on selected diagonals. These algo- 
rithms work under the more general Assumption 14.11 which is however not equally easy 
to check. In Section using this block algorithm, we present some examples of kernels 
for which Assumption 12.11 fails but nevertheless a perfect simulation algorithm can still 
be constructed. 
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2 Coupling of unilateral fields 

In this section we introduce the auxiliary i.i.d. fields which allow to couple, i.e. to 
represent on the same probability space, the family of unilateral fields X x defined in (J2J), 
for all possible boundary conditions x E E N+ x E N+ . Moreover we investigate when the 
values of these fields indicate that the dependence from the boundary conditions of the 
field in some fixed box is lost. 

In the next two sections we suppose that the kernel K satisfies the following minoriza- 
tion condition. 

Assumption 2.1. There exists a probability measure on E and a positive constant 
5 > So such that 

K{A\ yi ,y 2 )>6<t>{A) (6) 

for every measurable set A and any pair yi,y2 E E. 5o is a positive constant which will 
be specified in the next section. 

If 5 = 1 the field is trivially i.i.d. so we exclude this possibility from now on. 
Remark 2.2. If E is finite or countable, then provided 

T {z)=inf{K({z}\ yi ,y 2 ): Vl ,y 2 E E} (7) 
is not identically zero, the minorization condition (G|) is satisfied with 

5 = J2r(z)>0, (8) 

and 

*({*}) = (9) 
Under Assumption 12. II it is immediately checked that 

K{-\ yi ,y 2 )-6<t>(-) 



12/1,2/2) 
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is a kernel on E. Since E is Borel we can always define a function / : (0, 1) x E 2 which is 
separately measurable in each of its two arguments (see P) with the property that when 
U is uniformly distributed in the interval (0, 1), then f(U;yi,y 2 ) has the law H(-\y 1 ,y 2 ) 
for any pair (yi,y 2 ) € E 2 . Next we can prove the following 

Lemma 2.3. Let {Z, V, U} be mutually independent random variables with laws 

P(Z = 0) = 1-P(Z = 1) =5, V~(f), U uniform on (0,1), (10) 

and define 

g(z, v, u; y h y 2 ) = (1 - z)v + zf(u; y u y 2 ). (11) 

Then for any y±,y 2 E E 

g(Z,V,U; yi ,y 2 ) ~ K{-\ yi ,y 2 ). (12) 
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Proof. For any yi, y 2 G E we can write the kernel as a mixture 

K(-\ Vl , y 2 ) = 6<f>(-) + (1 - 6)H{-\ yi , y 2 ) (13) 

which is seen to be induced by the application of (fTTf to {Z, V, U}. □ 

We say that the family of functions {#(•; yi, y 2 ), (2/1, 2/2) G -E 2 } realizes a coupling 
of the family of laws 2/2) , (1/1,1/2) £ -K 2 } on the probability space where Z, V, 

and U are defined. Notice that when Z = all the random variables g(Z,V,U;yi,y 2 ) 
take the same value V^, irrespectively of (yi,y 2 ) G -E 2 - In this case we say that coupling 
occurs. When E is finite or countable, with 5 and <fi defined as in Remark 12. 2\ the above 
coupling is maximal, since the probability of coupling can never exceed 5. In fact, if 
{X yuy \ (2/1, 2/2) G -E 2 } is any other coupling of {K(-\y u y 2 ), {yi,y 2 ) G E 2 }, then 

P(3x g £ : X^ 2 = x, V(yi,y 2 ) G £ 2 ) < V inf P(X yi ' y2 =x) = 5. (14) 

Now let 

{Z, U, V} = {Zij, Uij, ^j}(ij)ez 2 
be mutually independent i.i.d. fields with (Zj j,Ui j,Vi j) ~ (Z,U,V) distributed as in 

From the previous lemma we see that the probability space supporting {Z, U, V} 
allows a coupling of the fields X x for all boundary conditions x G E N+ x E N+ by using 
the recursion 

Xfj = y(Z ; . r Vij, C,.f. X? .V, x ; ; ). i, j = 1, 2, . . . (15) 
along any increasing order of N +2 starting from the boundary conditions 

Xf Q = x ifl , = x 0ti i = 1, 2, — (16) 

The following definition will be useful in the sequel. 

Definition 2.4. For any subset S C Z 2 z£s external boundary, indicated with d S, is 
the set of parents of some elements of S, which are not themselves in S. Finally define 
the internal boundary of S , indicated with djS, as the set of vertices in S with at least 
one parent in d S . 

Notice that the external boundary of the finite box A m ^ n is the set of sites {(i, 0), (0, = 
1, . . . ,m, j = l,...,n}. 

It is immediately seen that we can use the representation (JT5jl to couple the field in 
any finite subset S of Z 2 for all possible boundary conditions prescribed on its external 
boundary d S, by recursion along any increasing order of the sites in S. For any choice of 
boundary conditions x = {xij G E, G d S} the unilateral field X x on S, with b.c.'s 
x and transition kernel K, can be therefore represented as 

X£ = G s (Z 5 ,V 5 ,U 5 ;x) (17) 
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where the function G s is suitably defined. 

Furthermore we denote by G S B {Z S) XJ S , V 5 ; x) the projection of the vector-valued func- 
tion G s (Zs, Us, Vs; x) on the sites belonging to the subset B C S. 

By the direct inspection of (fTTj) we notice that if Zq z s = then the recursion of (fTKj) 

can be started with the configuration Vg /( g irrespectively of the b.c.'s x G E as . We 
denote by V s {Zs, Us, V5) the resulting configuration of the field on S: hence 

Z dlS = => G S (Z S , U s , \ s ; x) = V s (Z s , Us, V s ). (18) 

The next lemma generalizes this situation. 

Lemma 2.5. Lei S be a finite subset of 1? , and B C S . IfZg^ = 0, then 

G s B (Z s ,U s ,V s ;x) = r B (Z B ,U B ,V B ), Vxe-B^. (19) 

Proof. By the argument above the restriction of the field on B does not depend on its 
values on d B. By consequence it does not depend on the b.c.'s x on d S. □ 

3 Oriented percolation and coupling 

In this section the coupling made above is exploited to construct, under Assumption 12.11 
a perfect simulation algorithm for the equilibrium distribution of an ergodic unilateral 
field in any finite box A m>n . Since we keep the box fixed, we write it simply as A. Instead, 
the boundary conditions are pushed far away from A, and the field is constructed using 
a finite number of samples of a given realization of the auxiliary fields Z, V and U. 
In particular the auxiliary field Z will play a crucial role, since it indicates when the 
construction becomes insensitive to the boundary conditions. 

We begin by defining an (increasing) path 7 joining two vertices as a sequence of 
vertices {(ik,jk) £ 2?}fc=o,i,...,m> with m > 0, such that for k = 0, . . . , m — 1 either 

ik+i = ik + l,3k+i=3k 01 i fe+ i = ik,jk+i = 3k + 1- (20) 

In this case we say that (i , j ) and (i m , j m ) are joined by the path 7 of length m. Obviously 
two distinct vertices (a, b) G Z 2 and (c, d) G Z 2 can be joined by a path if and only if 
a < c and b < d. 

Given a realization of the field Z = {Z it j G {0, 1}, G Z 2 } we say that a path 
7 = {(zfe, jk)}k=o,i,...,m is open (in the field Z) if Z ikt j k — 1 for k — 1, . . . , m. By convention 
we do not require that Z io> j = 1. 

We are now interested to study the following random subset of Z 2 

00(A) = u(A, Z) = G Z 2 : 3 an open path in Z joining to (k, I) G <9/A} = C k ,i 

(fc,Z)69/A 

(21) 
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Figure 1: a particular realization of a; (A). 



where Ck,i is the set of vertices joined to (k, I) G dj A by an open path in Z (see Figure [T] 
for a particular realization). 

A fundamental result of oriented percolation for Bernoulli fields (for a general reference 
on this subject see pQ) ensures that if the probability of 1 at a site does not exceed a 
critical value < p c < 1 then j is finite with probability 1 (irrespectively of (k, I) 
by translation invariance). Moreover if such probability is strictly smaller than p c then 
E|C^| = E|C 0j o| < oo (see again [T]), which implies that 

E|t^(A)| < |a/A| -E|C 0j o|- (22) 
Since an estimate from below is trivially obtained using the bound 

KA)|> Yl z ^ ( 23 ) 

(M)e£>/A 

we conclude in this case that E|u;(A)| = 0(|<9/A|). 

From now on we set 5q = 1 — p c in Assumption 12.11 Thus 

Proposition 3.1. Under Assumption \2. II the random subset u(A) is finite almost surely. 
Moreover if 5 > 5q then the mean value of \uj(A)\ is 0(\djA\). 

Since we are interested in bounding 5 from above, we need to bound p c from below. 
For example in [7] it is proved that p c > 0.682 . . . , thus for 5 > 0.317 . . . Assumption 12.11 
is verified. 

Let us define the random set B(A) = A U u(A). Then we can prove the following 
Lemma 3.2. For any {i,j) G diB(A) it is Zij = 0. 
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Proof. We first notice that the internal boundary djB(A) is the disjoint union of '■ 
{i,j) G <9/A, Z it j = 0} and diu(A). Therefore, to prove the lemma, we need only to show 
that Zij = for any (i, j) G diuj(A). This is due to the fact that if Zij = 1 then both the 
parents (i— 1, j) and (i, j — 1) are in u(A) since the existing open path from (i, j) to a site 
in diA can be extended in both directions. But this is absurd, since (i, j) G diu(A). □ 

At this point the main result of the paper can be proved. As before the function 
denotes the restriction of T B to the sites in A C B. 

Theorem 3.3. Under Assumvtion \2.1\ the kernel K is uniformly ergodic. Moreover for 
any finite box A = A m .„ 

rf (A) (Z B(A) , V fl(A) ,U B (A)) ~ X' A , (24) 
where X' is the equilibrium field with kernel K . 

Proof. By Proposition 13.11 B = B(A) is almost surely finite. Next, for any finite S such 
that Ac S, let Q$ be the event {B(A) C S}. From Lemma \2. 51 and Lemma 13. 21 we know 

that if Qg occurs, for any x 6 E 9S 

Xg |A = G S A (Z S , V 5 , U 5 ; x) = rf (Z fl , Vs, U B ). (25) 

Next, let {Sk D A,k = 1,2,...} be any sequence of finite subsets increasing to the 
set of all sites (h, I) G Z 2 such that h < i and j < I. Then, under Assumption 12.11 the 
sequence {Qs k } increases to an event of probability 1. This means that with probability 
1 B(A) C Sk eventually with k, hence 

G s A k (Z Sk , V Sfe , U Sk ; x fc ) = rf (Z B , V B> U B ) Vx fc G E^ 5fc (26) 

eventually in k. 

Next let us consider any sequence {(hk,Vk)k = 1,...} converging towards oo and 
construct the sequence of boxes 

C k := {1, . . • , h k + m} x {1, . . . , v k + n} C (/i fc , u*,) + A =: A fc . 

It is clear that for any b.c. x^ G E 9Ck 

X (H k ,v k)+ A = G C A :(Z Ck , V Cfc , \J Ck - xfc) ~ Gt(Z Sfe , V Sfe , U Sfe ; x fe ) (27) 

where 

Sfc = {-/i* + 1, • • • , m} x {-w fc + 1, . . . , n}, 
with the b.c.'s at the r.h.s. of ()27|) defined on 

9 S k = {{-h k , -v k + 1), . . . , (-h k , m)} [J{(-h k + 1, -v k ), . . . , (m, -v k )}. 
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The r.h.s of (|27j) is equal to T^(Zb, Vb,Ub) on Qg . By consequence, denoting by 1a 
the indicator of A 

IIA^ fc)+A )-Arf)il = ||£(G?)-£(rf)|| = supE(|/ A (G?*)-/ A (rf)|) < 2P(Q£/), 

(28) 

from which we deduce uniform convergence in variation because lim^oo P{Qg, ) = 0. 
Since, by varying A, this specifies a stationary and compatible family, by Kolmogorov's 
Theorem it is a realization of a stationary field X' in the finite box A. □ 

To summarize we have constructed a sampling scheme for the equilibrium field X' 
on a finite box A, based on a random but a.s. finite number of samples from the fields 
Z, V, U. In fact, by Proposition 13.11 this number is proportional to |A| plus 0(|<9/A|) in 
the average. 

We conclude the section by discussing a possible implementation of the algorithm. We 
construct uj(A) by backward induction in the following way. Let 

A = {(i,j):(i,j)ed I A,Z hJ = l} 

Then for k > 1 we determine A&+i from A& as 

A fe +i = : {(i + + 1)} D A k ± 0, Z id = 1} 

until the index k max such that Ak max +i = for the first time, which is finite with prob- 
ability 1 if 5 < 5q- Then oj(A) is the union of U^^A^ and its external boundary. Then 
order the sites in a; (A) by starting with those in such an external boundary, then those in 
Ak max , next those in Ak max -\ which were not already in Ak max , and so on. After having 
totally ordered all sites in a; (A), we start with those in A. No matter which order is chosen 
within each of these regions, the recursion along this total order will allow to compute 

rf(z B ,v B ,u B ). 

4 A class of more general block algorithms 

In the previous section we have shown how to construct a perfect simulation scheme for 
some stationary ergodic unilateral fields on Z 2 . We have used percolation arguments on a 
suitable auxiliary Bernoulli field to show that the algorithm works under the " sufficiently 
large" minorization condition Assumption 12. II On the other hand, in the one-dimensional 
case the same idea leads to an algorithm for the exact simulation of a sample from a 
stationary discrete-time Markov chain that works under any non trivial minorization 
condition, since in any non degenerate one-dimensional Bernoulli field clusters are always 
finite. 

More generally, the Multigamma coupler of Murdoch and Green extends the above 
idea to cover the whole class of uniformly ergodic kernels by considering a suitable power of 
the kernel. In fact, by Theorem 16.0.2 in jHj uniformly ergodic kernels are characterized 
by a minorization condition on some power K m of the kernel K, which means that we can 
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apply the same algorithm to the m-skeleton chain {Xk m , k = 0, 1, . . .}. Moreover provided 
m is large enough the value of 5 in the minorization condition can be taken arbitrarily 
close to 1. More generally, it has been proved in [Sj that a vertical backward coupling 
time exists only for such a class of Markov chains. 

Based on these observations, in this section we try to extend the previous results to 
a wider class of unilateral fields by considering suitable skeleton fields, which in general 
consist of blocks of sites. 

By a diagonal of Z 2 we mean a set 

D h = {(i,j)eZ 2 :i+j = h} 

for h G Z. The distance between Dh and Dk is defined as \h — k\. The binary field used to 
indicate the region where the simulation have to be performed is in general defined over a 
new lattice associated to blocks of I adjacent sites lying on diagonals at distance (d — 1)1 
one from the other, for given integers 1 = 1,2,... and d = 2,3, . . .. We start by defining 
the block 

B 0>0 = {(1,1), (2,1-1)..., (1,1)} (29) 
and, for any pair (i, j) G Z 2 such that i + j = mod(d — 1), define its translates 



Bij = B 0fi + (il,jl) 



(30) 



B 



B, 



B, 



u-i 















B l,-2 





R 



i.-i 



Figure 2: Blocks with / = 2 and d = 2. 



These blocks are taken to be the vertices of a new graph Q d = (V^,^). An oriented 
edge connects Bh,k with Bh+i,k+d~i-i for i = 0, 1, . . . , d — 1, for any pair (h, k) G Z 2 such 
that h + k = mod(d — 1). In the following we will directly identify Vd with such a subset 
of sites of Z 2 and refer to the sites (h + i, k — i — d + 1) G Vd as the d parents of the site 
(h, k) G V d , for i = 0, 1, . . . , d — 1. It is clear that for I = 1, d = 2 we get that the graph 
Qd is the original lattice Z 2 ; more generally, for larger values of I and d = 2, the graph Qd 
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is isomorphic to Z 2 . In the next section we will devote Example 2 to show that it may be 
convenient to let I grow, rather than d. 

Next let S = {s m ,m = 0, 1, . . ., (d — 1)1} be a family of coupling functions, i.e. for 
m = 0, 1, . . . , (d — 1)1, s m : (0, 1) x E 2 — > E is separately measurable in each of its 
arguments and such that 

s m (U;x 1 ,x 2 ) ~ K(-\x 1 ,x 2 ), (31) 

when U is a uniformly distributed random variable taking values in (0,1). Since E is 
Borel the function g defined in flllj) is a particular example of ()31|). since all the auxiliary 
random variables Z, V and U which are required can be constructed as functions of a 
single uniformly distributed random variable. The construction of the field X* with kernel 
K over any finite region S with b.c. x on E 9 s is then performed by using a vector Us of 
independent uniformly distributed random variables, iterating along an increasing total 
order of the sites in S the recursion 

Xfj = s m (i,j) (U itj ; X*_ l p Xf^) (32) 

where 

m(i,j) = % +j — I — 1, mod{(d—l)l}, 

starting from 

X h = x - 

Notice that the coupling function s m used to construct the value of the field at a given site 
is allowed to depend on m(i,j), the distance of the diagonal where lies from 
the "previous" diagonal of blocks. In Example 1 we will show the usefulness of allowing 
couplings depending on the diagonal. 

Next, for any fixed pair of integers (h, k) G V<f, consider the trapezoidal region 

Rh,k = {(ij) ■ i < lih + < l(k + 1), (h + k - d + 2)1 + 2 < % + j < (h + k + 1)/ + 1}. 

The reason for defining this region is that Bh,k C Rh,k and 

d Rh,k = ^i=oBh~i,k+i-d+i- 

Therefore we may represent the field as a function of the vector Un hk with i.i.d. 

components and the boundary values on the parent blocks 

x= (x Bh _ ik+i _ d+1 ,i = 0,...,d-l) 

in particular 

X^ fe = F 5 (U Rhifc ;x). (33) 

where F s is defined through the recursive application of (j32j) along an increasing total 
order of the sites in Rh,k- 

We are now ready to make the following general assumption. 
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Assumption 4.1. There exists a binary field 

W h , k =^(U RhM ), (h,k)EV d 

with 

P(W h , k = 0) = ~5> ±Z±, (34) 

a 

such that {Wh t k = 0} implies 

F s (V Rh!k ; x) = $(U Rh ,J, Vx G E® Rh ' k (35) 
for some measurable function $. 

We kept Assumption 14.11 quite general in order to accomodate various possible defini- 
tion of the field W = {W^k, {h, k) G V^}, for a given family S of coupling functions. 
In principle (at least when the state space E is countable) the field W can be directly 

defined to have the value zero if and only if F s (Un h k ; x) does not depend on x G E dRh - k , 
and in this case §(XJf> hk ) is equal to such a common value. This means that in order to 
conclude that Wh,k = we have to check that the realizations of the field k started 

from all the possible b.c.'s x G E dRh > k collapse into a single value. Thus, within a 
single region Rh, k this is similar to the original Propp and Wilson coupling from the past 
algorithm [TH] . 

A computationally less demanding choice is to define recursively the random subset 
of E 

= Sm(i,j)(Uij', Iij-i X Jj_]_ ( j-) ; (36) 

for any site in (i,j) G Rh,k, starting from Jjj = E for G d Rh,k- Finally we define 
Wh,k = if Iij is a singleton for all G Bh,k, and $(Ui? hfc ) is then equal to its unique 
element. Another option for the definition of W will be presented in Example 1. 

A generalization of the coupling is possible by allowing Sm(ij)(s -^f-iji X? to 
depend on Ui >n , (I, n) G Rh,k, with I < % and n < j, in such a way that 

s m (i,j){Ui >n , (l,n) G R h , k ,l< i,n< j; ~ K(-\X*_ l;j , X*^). 

A particular example is when s m ^j) realizes the maximal coupling of the laws {K(-\x\, X2), 
(xx,X2) G x Iij-i} where I^j is defined as in (JBTIj) . 

Now consider any finite subset 5 C Z 2 . By analogy with Lemma 12.51 it is quite clear 
that if the set B C S is such that diB is contained in a union of blocks B^^ where 
Wh^k = 0, then the restriction of the field X£ on B can be represented as a measurable 

function Y B (XJ B ), irrespectively of the b.c.'s x G d S. 

We say that 7 = {(i k ,jk) £ Vd}fc=o,i,-,ro is an increasing path of length m if for 
k — 0, . . . ,m — 1 the vertex (ik,jk) is a parent of jfc+i) in the graph The field 
W is then used to define open paths in 
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Notice that a field Wh,k — ^(U.R hfc ) defined for (h,k) G Vd (see Assumption 14. 1ft is 
not a Bernoulli field in general. However since U is an i.i.d. field the random variables 

{W hiM ,(hi,ki) G Vd, i = l,...,t} 

are mutually independent whenever the regions Rh i7 ki are pairwise disjoint, for 2 = 1,...,*. 
Now it can be verified that RhiM an d Rh 2 M are n °t disjoint if and only if h\ + k\ = 
h 2 + k 2 = j(d — 1) for some integer j and \hi — h 2 \ < d — 1 (e.g. check in Figure El that 
Rofif)Ri-i = {(2,0)}). Therefore the field on each diagonal is (d — Independent and 
different diagonals are mutually independent. 

For the sake of simplicity, in this section we suppose that djA is a union of blocks 
Bh t k, so it can be identified with a subset of Vd] otherwise we can enlarge A to have this 
property. In analogy with (j21j) we define the random region 

cr(A) = <t(A, W) = G Vd '■ 3 an open path in Qd joining to (k, I) G <9/A}. 

(37) 

and D(A) = A U cr(A). As in Lemma f3. 21 we can prove that for any (h, k) G diD(A) the 
random variable Wh,k — 0. We can finally show the following 

Theorem 4.2. Suppose the kernel K satisfies Assumption Then K is uniformly 

ergodic and for any finite region A the law of its stationary version X'a is the same as 
T D ( A )(U D(A) ). Moreover E|Z?(A)| = |A| + 0(|9/A|). 

Proof. As in flU} 

^(A)= |J C' h>k 

(h,k)ediA 

where C' h k is the set of blocks which can be joined to (h, k) by an open path in Qd- We 
have only to prove that under Assumption 14.11 C' Q is finite almost surely and has finite 
mean. The only difference with Theorem 13.31 is that W is not Bernoulli. However, since 
different diagonals are independent, the values of the field W on any increasing path are 
i.i.d. Next define 

diam(Coo) = max{ length of an open path from (h, k) to (0, 0): (h, k) G C 00 }. (38) 

Then the following holds 

P(K >0 \ > k 2 ) < P(diam(^ >0 ) > k) < {d(l - 5)}\ (39) 

where the first inequality is due to a trivial geometric argument. The second inequality 
follows since the probability that a fixed increasing path of length k is open is equal to 
(1 — S) k and the number of such paths leading to the origin is equal to d k . Since by 
Assumption 14.11 d(l — 5) < 1 we obtain that C' 00 is finite w.p. 1 and its cardinality has 
finite mean (in fact, finite moments of any order). □ 
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It is natural to ask whether it is possible to improve the inequality (}3"4*j) by dominating 
W with a suitable Bernoulli field. A domination result of this type for (d — Independent 
fields can be found in [TT] : however we have checked that it does not improve the bound 

dSH). 

Before discussing some examples we need to make clear the limitations of the block 
algorithm presented in this section, due to the difficulties arising from the two-dimensional 
structure. In some sense we have constructed a skeleton process , (h, k) G Vd' 

however, this field is not a ^-unilateral field anymore, due to the overlap between the 
regions Rh,k, for (h, k) e Vd- As a consequence a minorization condition of the type 

P(X£ G A\X^ = x) > S(p(A) 

y B h,k I dR hik 

for all x e E dRh < k , does not immediately translates in the possibility of coupling the field 
with an auxiliary Bernoulli field with probability 5 of 0. For this reason we need the 
stronger Assumption 14. 11 which requires the choice of a specific family of functions {s m } 
allowing to simulate the field also at sites common to different regions. Moreover, since 
the binary auxiliary field used by the algorithm is not Bernoulli, the values of 8 for which 
we can prove that the algorithm works have to be larger than those obtained with i.i.d. 
percolation (5 > 1/2 rather than the previously cited bound 5 > 0.317 . . . , for d = 2). 

As a consequence we cannot identify the class of uniformly ergodic fields which satisfy 
Assumption 14.11 for some choice of / and d. As a matter of fact we are not aware of 
any simple characterizations of the class of uniformly ergodic unilateral random fields, as 
possible in the one-dimensional case. 

5 Examples 

Example 1. In this example we will construct a perfect simulation algorithm for the 
stationary unilateral field with kernel K in a finite region under the following 

Assumption 5.1. There exists a measurable subset C C E, a probability measure on 
E and positive constants < p\ < 1 and < p<i < 1 with p\p 2 > 2/3 such that 

K(C\y 1 ,y 2 )>p 1 W( yi ,y 2 )eE 2 (40) 

and 

K(A\ yi , y 2 ) > P 2<P(A) V (y h y 2 ) e C 2 (41) 
for every measurable set A. 

Under Assumption ^, ll we will build a coupling of the field with kernel K which satisfies 
Assumption 14.11 with I = 1 and d = 3. In this particular case we need to distinguish 
between even and odd numbered diagonals, depending on the parity of the sum of the 
coordinates of a site: the blocks will be single vertices (h, k) lying on even numbered 
diagonals. 
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The field will be constructed by means of two families of functions so and s\ with the 
property (pTTj). which are used to get the value at sites lying on even and odd numbered 
diagonals, respectively. For simplicity of notation we use pairs of random variables (U\, U 2 ) 
which are uniformly distributed in (0, l) 2 rather than a single random variable in (0, 1). 
The functions sq and s\ are defined by means of some functions /$ : (0, 1) x E 2 — > E, for 
i — 0, 1,2,3,4, measurable in each of the arguments, with the following properties. For 
(2/1,2/2) g£ 2 

fo(U 1 ',y 1 ,y 2 )=f (U 1 )^ < f ) (.), (42) 

fl{u . M2) „ J i^TOifc.i*) - **(•)] for m,y 2 g C 

I K (-|j/i, y 2 ) otherwise. 

h{U X ] yi, y 2 ) ~K{- \ yi ,y 2 ) . (44) 

f (tt \ K i.-^C\y 1) y 2 ) . . 

fsiUiWuyz) 7T 77 r l v— , 45 

K{C\y 1 ,y 2 ) 

f itt x 1 x K(-nC\ yi ,y 2 ) \ 

h{Ui)yi,y 2 ) ~ <^ K(-\ yi ,y 2 - p x — \ . 46 

1 - Pi I K(C\y u y 2 ) J 

These functions always exist since E is assumed to be Borel. Notice that (|46|) is well 
defined since 

myi,y2)-Pi K ^J: lm,y , 2) = K(.nc c \ yi ,y 2 )+(i - — -) K(.nC\ yi ,y 2 ) (47) 

^(C lj/1,2/2) V K[C 12/1,2/2)/ 

is non-negative by (pfUj) . 

With these positions it is not difficult to verify that the functions 

„ („ m • ni it \ \ ( X ~ 1 (o,p 3 )( u 2))/o(wi) + 1(o, P2 )(m2)/i(mi;2/i,2/2) for 2/i,2/ 2 G C, 

*o(«i,«2,Vi,lfe)-| /2(tti;yijlte) otherwise. 

(48) 

Si(z,Ui,U2;2/i,y 2 ) = 1 (o,pi)( w 2)/3(wi; 2/i, 2/2) + (1 - l(o,pi)(w2))/4(wi; 2/i, 2/2), (49) 

satisfy t)31j). Essentially so and si come from two different mixture decompositions of the 
kernel. Finally the field W is defined for any pair (h, k) G Z 2 such that (h + k) is even 
by 

W h , k = 1 " l(0,p 1 )(( f/ 2) /l -l, fc )l(0,p 1 )((^2)M-l)l(0,p 2 )((^2)M)- (50) 

Notice that Wh,k = if and only if all the indicators appearing in the above formula are 
equal to 1, which implies that Xh,k = fo((Ui)h,k), irrespectively of the values at the parent 
sites. Since 

P(Wi,j = 0) = P \ P2 > 2 - 

we have thus established 
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Theorem 5.2. Under Assumption \5.1\ the conclusions of Theorem \4-ty hold. 



It is not difficult to construct a kernel K for which Assumption 15 . II is satisfied, whereas 
Assumption 12. II is false. Let E = {0, 1, 2} and suppose that 



Choosing C = {0, 1} gQ| holds with Pl = 0(C) = 0(0) +0(1) and (jUJ) holds with p 2 = 1, 
hence if 0(C) is large enough Assumption 15. II is satisfied. On the other hand Assumption 
12. II does not hold, since K(-\2,2) and i^(-|l,2) are singular. 

Example 2. This example allows a simple analysis of the role of d and I in the block 
algorithm. On the state space E = {0, 1, 2} consider the unilateral kernel 



It is clear that this kernel does not satisfy Assumption 12.11 for any p 6 [0,1], since 
K(0\2,2) = #(1|1,1) = AT(2|0,0) = 0. Now assume p < 1/2 and consider what hap- 
pens by increasing d or I. For d = 3 and I = 1 we see that by choosing C = {0, 1} and 
the Dirac mass on 0, with pi = p 2 = 1 — p, we can fulfil Assumption 15.11 provided 
(1 — p) 3 > 2/3. This choice of C is clearly the best possible. 

Next we show that by taking I = 2 and d = 2 we can enlarge the region of parameters 
p under which uniform ergodicity can be proved and a perfect simulation algorithm can 
be constructed through Assumption 14.11 For any site of the lattice we always use the 
natural coupling s : (0, 1) x E 2 — > E, defined by 

s(u;y 1 ,y 2 ) = l(o, P )(«) min{2, min(yi, y 2 ) + 1} + l(p,i)(u) max{0, min^, y 2 ) - 1}. (53) 

We define W ,o = provided £/y > p on all the sites = (1, 2), (1, 1), (2, 1). The whole 
field {W/i,*;, (h, k) G Z 2 } is defined by translation. Such a definition ensures that the field 
has the value zero on both the sites (2, 1) and (1, 2) which form the block -Bo,o, irrespec- 
tively of the b.c.'s on the parent blocks, formed by the sites (—1, 2), (0, 1), (1, 0), (2, — 1) 
(see again Figure |2J). Hence we can take $ = in (|35|). It is immediately obtained that 
P(W fl = 0) = (1 — p) 3 , so that Assumption 14.11 is satisfied provided (1 — p) 3 > 1/2, 
enlarging the region where we it is proved that perfect simulation works. 
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